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ABSTRACT 


Several experimental and theoretical studies have shown that path>integrated rain rates can 
be deic nnined by means of a direct measurement of attenuation. For ground based radars this is 
done by measuring the backscatteiing cross section of a fixed target in tire priisence and absence 
of rain along the radar beam. A ratio of the two measurements yields a factor proportional to the 
attenuation from which the average rain rate can be deduced. In tfiis paper, we extend the tech' 
nique to spacebome radars by choosing the ground as reference target. The technique is also gen* 
eralized so that both the average and range-profiled rain rates can be determined. The accuracies of 
the resulting estimates are evaluated for a narrow beam radar located on a low earth orbiting 


satellite. 
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DETERMINATION OF RAIN RATE FROM A SPACEBORNE RADAR USING 
MEASUREMENTS OF TOTAL ATTENUATION 

INTRODUCTION 

Radar is the most promisiiig spacebome instnunent for quantitative measurements of rainfall 
rate over land. For this q>plication, however, radars that operate at non*attenuating wavelengths 
have several disadvantages. The restrictions on antenna size from space preclude fields of view 
small enough at the longer wavelengths to resolve the more intense regions of rain rate. As the 
resolution requirements are relaxed, the partially filled range cells lead to errors caused by spatial 
averaging, small signal to noise ratios, and reflectivity gradients. Apart from the errors caused by 
poor resolution, the Z-R law (used in deriving the rain rate R from the reflectivity factor Z) is a 
sensitive function of the unknown drop size distribution (DSD) [ 1 ] . When conditioned on a 
knowledge of the synoptic storm conditions, some of the large scale variability in the DSD can be 
eliminated; nevertheless, a significant variability persists within each class of storm [ 2] . The result- 
ing errors, along with tyjiical calibration errors of 2db to 3db limit the accuracy of rain rate deter- 
mination from a single -iOn-attenuating wavelength radar. 

While short wavelength attenuating radars offer better resolution for a fixed antenna size, the 
measurement accuracy of R is degraded not only by errors in the calibration constant and the Z-R 
relationship but also by errors in the attenuation coefficient-rain rate law (k-R). As in the Z-R law, 
a major source of error in the k-R law arises from systematic changes in the DSD as well as fluctua- 
tions in the DSD within each storm class. For moderate to large values of attenuation, even small 
errors in the k-R relationship can introduce a large bias into the rain rate estimate [3] . 

The addition of an independent measurement of some meteorological quantity can often 
reduce one or more of the errors mentioned above. One of the first examples of this type of correc- 
tion procedure was given by Hitschfeld and Bordan [31 . These authors proposed an independent 
measurement of rain rate at some point along the radar beam which served to bound the errors in 
the rain rate estimate. The dual-wavelength methods [4-6] can also be viewed as a type of correc- 
tion procedure in which the measurements arc used to derive a two parameter DSD. The liquid 
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water content and the rain rate can then be deduced without the need of Z*R or k*R laws. Other 
dual measurement techniques include the use of orthogonal polarizations [7,8] and the use of a 
combined radar-radiometric sensor [9]. 

Of particular importance for the work presented here is the method [ 10-14] in which a direct 
measurement of total attenuation is used to find the path-integrated or the path-averaged rain rate. 
The attenuation is obtained from a ratio of iJie backscattered powers from a fixed target where one 
measurement is made in the presence, the other in the absence of precipitation. From the attenua- 
tion, an appropriate k-R law, anu t>^ ^ path length from the radar to the target, the path-integrated 
or average rain rate can be detemiip.;id without recourse to the Z-R law or the radar calibration 
constant. Atlas and Ulbrich [ 14] have shown, moreover, that at w&velengths near 0.86 cm., tfie k-R 
law is essentially independent of the DSD. This independence is manifested in less scatter about 
the k-R regression line and therefore a reduction in the error variance of R. 

For satellite applications, the technique is not directly applicable because of the impractic- 
ability of a global network of calibrated targets. The objectives of this paper are twofold: to modify 
the technique to make it suitable for spacebome radars and to generalize it so that, in addition to 
the average rain rate, range-profiled rain rates can be estimated. 



THE SURFACE REFERENCE TECHNIQUE 

By means of an attenuating-wavelength spacebome radar, measurements of the backscatteiing 
coefficient (i.e. the scattering cross section per unit area) of the surface, o° , can be made in the 
rain-free areas adjacent to the rain volume. As the beam passes into regions of precipitation, the 
backscattered power from the surface will be proportional to o^A where A is the attenuation 
factor. If the incidence angle and the radar frequency can be chosen so that the values of a** within 
and outside the precipitation are nearty invariant, then the ratio of the two measurements yields an 
accurate estimate of A. Once A has been obtained the average rain rate follows from a k-R law and 
a measurement of the path length through the rain [IS]. 

In the s:andard methods, is obtained from a measurement made in the absence of precipi- 
tation. In the present technique, which we will refer to as a surface reference technique (SRT), the 
actual 0 ° is an unknown which can be approximated either by an average value measured in the 
rain-free areas adjacent to the target or by a subsequent measurement of at the actual target 
when precipitation along the radar beam is absent. The replacement of the actual with an 'equival- 
ent* reference measurement introduces an additional error into the determination of rain rate, the 
magnitude of which depends on the spatial or temporal inhomogeneity of the surface scattering 


properties. 

To restate the above ideas in a more quantitative maimer, we let denote the radar return 
power from the rain at the jth range bin. We further classify the surface return according to 
whether precipitation is present or absent along the radar beam, choosing as the surface return 

in the presence of precipitation and P^. as the unattenuated surface return. In the subsequent 

j 

discussions we will often have need to distinguish between those range bins (within the main beam) 
which intercept the surface from those which do not. We refer to the former as surface bins and the 
latter as precipitation bins as shown in Fig. 1. 

As shown in the appendix the return powers at the nth range bin can be approximated by 
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( 2 ) 


® 2 * 

where 

f*“ 

-0.46 j k dt 

\aje '® (3) 

k is the attenuation coefficient in db/km, is a calibration constant, r^ is the distance from the 

radar to the edge of the storm, and r^ is the range corresponding to the center of the nth range bin. 

The subscripts on are used to indicate that the values of o'* for the two measurements are 

generally different The ratio of P'q toP^ then provides an estimate of A„. In general, however, 

n n 

the return power from the surface cell, n, in the presence of precipitation will include a contribu> 
tion from this precipitation (see Fig. 1). If we also account for the noise power, , added to the 
return by the receiver, then the estimate of written A„ , is given by the ratio of the total return 
powers measured at the surface bin n in the presence and absence of precipitation: 

( 4 ) 


where 






(5) 

( 6 ) 


The subscripts added to P^^ indicate that these contributions are not constant over the two measure- 
ments. 

Using the k-R law, k» 7 R^, in (3) and the estimated value P^ n^ 2 n 




( 7 ) 


Defining the average rain rate to be 


1 f “ 


( 8 ) 
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with 


r,-t.-t, (9) 

we notice diet for an arbitrary spatial distribution of rain rate along the beam and | # 1 , can 

not be expressed in temis of Pjg, P 2 n. Instead, we approximate the average rain rate by 



1/1 


so that upon using (7), the estimated rain rate can be written in terms of the measured quantities 
Pj a* ^ 2 a> assumed values ( by 

r I T'* 

For wavelengths between 0.7 cm and 2 cm we can use the appro cimation | ■ 1 and adjust y 
accordingly without a great deal of increase in the variance in the k*R law [ 14] . 
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RANGE PROFILING OF THE RAIN RATE 


To obtain profiled rain rates (i.e. rain rates at each range bin) we can begin either with the 
Hitschfeld'Bordan algorithm [3] or the iterative algorithms proposed by Oimsby [ 16] . Although 
these techniqi’ .is themselves provide range profiled rain rates, they are susceptible to large errors 
whenever the total attenuation is significant. Since the SRT is feasible only when moderate to large 
attenuations occtir, the above methods are of limited applicability in this instance. We can, how- 
ever, reformulate them by using the total attenuation, obtained from the SRT, to bound the errors 
in the rain rate. In fact, the procedure to be described is similar to the Hitschfeld-Bordan correction 
technique mentioned above except that the attenuation estimate of (4) now takes the place of the 
independent measurement of rain rate. 

For the particular application used here, the Hitschfelu-Bcrdan algorithm has been shown to 
yield better results than the iterative estimates. We therefore deal exclusively with the former. 

The Hitschfeld-Bordan algorithm for the rainfall rate at the jth range bin can be written [ 17] , 


Rj«-» 


«aZ 


m, 



( 11 ) 


with 



i^j 

i“j 


The reflectivity factor Zj, the measured reflectivity factor Z^^ and the calibration constant 
related by the equations 


ate 




where r, P 




. A, have been defined previously. 


- 0.46 


Aj »e 


Using the k-Z law, k * oZ^ and (3) then 

at J ^ Z^it) di 


( 12 ) 

(13) 


(14) 
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Equation (11) follows from the Z-R relationship, R ■ aZ^ , equation ( 1 2) and the fact that 



which is a consequence of ( 1 2) and ( 14). The quantity h in ( 1 1 ) and ( 1 5) is the radar range lesolu* 
tion which arises in ^proximating an integral by a summation, i.e., 

|^2^(s)d8'*h 2 ■ 

We now assume that an estimate of the attenuation factor up to the nth range bin, A^, is available. 
There are several ways of incorporating this additional measurement into the rain rate estimate of 
(11). We discuss two possibilities. 

Case 1 : To account for an offset in the calibration constant Cj^ , we multiply the set |z^ 
i» 1,2, ... , jin (11) by some constant, say p. Since Z^^ » r^^/C^, this operation is equivalent 

to multiplying C,^ by p'^ . Equation (11) then becomes 

( J 

l-0.46o^h £ (16) 

Identifying pZ^ ^ as a corrected form of Z^ then from i 1 S), 

[ i 

l-0.46«Sh 2 «,(pZ„^)?J _ 

A 

Letting j increase to n in this formula, using for the estimated value A^ from (4), and solving 
for p givvs 

[ « g 1 

(l-A,®)/0.46a^h I «,Z"^ 

The quantity p is now inserted back into (16) which yields a new rain rate estimate that we label 
with a superscript 1 : 
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vheie 



(17) 

J • (( 1 - A,**) f 0.4« « # 

(18) 


(19) 


i*i 


These equations are valid for j ■ 1 , ...» n. Equation (17) can be interpreted as an alternative form of 
the Hitschfeld-fiordan algorithm in which the quantity has been used to eliminate offset errors 
in the radar calibration constant, C,^ . 

Case 2: To compensate for offset errors in ot, where k ■ azfi, we replace a by op' in ( 1 1) and 
(IS). Proceeding in the same way as in Case I, we find 


p'-(l-A„^/0.46a;?S„ 

The quantity p', in turn, yields a second rain rate estimate, that can be written 

(1 -(1 - A.") s, / (20) 

forj^ n. 

Although the only difference between ( 1 7) and (20) is the factor of J appearing in the former, 
it can be seen upon inspection that (20) is independent of offsets in a but depoident on errors in 
the calibration constant. Exactly the opposite is true of ( 1 7). 

In both cases considered the total attenuation, as determined by the SRT, was used to derive a 
parameter which provide; a bound to the Hitschfeld*Bordan algorithm thereby eliminating the 
occurance of extremely large errors. It is worthwhile noting that a number of other range profiled 
rain rate estimates can be deduced. For example, if the path integrated rain rate. R^. is defined to 
be R ds, then from (7) an estimate of this quantity is. 


■If 


T i/q 




( 21 ) 
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On the other hand, can be estimated from any of the following quantitim: 

h t *)"■■ . 

h t R,<‘> . 
h t . 

where ^ an given by (11), (17), and (20) respectively. 

If we multiply by some constant certain variables in any of these three quantities and adjust it until 
an equality with (21) is obtained, the constants found in this way supply new estimates of range pro* 
filed rain rate. In this paper only the estimates given by (10), (17), and (20) will be analyzed. An 
analysis of the other estimates would probably be warranted only after an experimental veriflcation 
of the basic technique. 
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RESULTS AND DISCUSSKW 

The remainder of the piQ)er will be spent in investigating the behavior of e<iubtions (10), (17) 
and (20) in the context of a down^ooking radar located on a low earth orbiting platform. Ln par* 
ticular, we compute the mean and standard deviation of y , \ ’ as functions of rain rate 

for suitable values of transmit power, wavelength, antenna diameter, etc. The parameten that will 
be kept constant ate listed in table 1. Since the quality of the \ R*^^ estimates varies markedly 
with distance into the storm, their statistics are presented as a Ainction of range. 

The procedure for finding the statistics of the above estimates can be conveniently described 
as a two>stage process. In the first stage, the average values of the return power from the rain and 
the surface are computed as a function of distance into the storm (see appendix). The second stage 
in this procedure is to introduce offset errors and random variations into several of the radar and 
meteorological quantities. Via a simulation, these quantities are used to compute particular realize 
tionsofft^y, andR^^^ inaccordantt with equations (10), (1 7), and (20). Grntinuing in this 
way over one thousand such realizations, we then compute the sample mean and standard deviation 
ofR^y andof R^‘^ R^^* forj* 1, ..., n where j ■ 1 corresponds to the first bin contaiiimg rain at 
the storm top while j « n corresponds to the bin at which the surface return power is maximum. 

An example of the mean power calculation is shown in figure 2 where the receiver noise 
power, Pj 4 (N), the rain return, Pj^ (R), and the surface or clutter return, P*q (C) have been plotted 
4t each range bin. The abscissa represents the slant range which is measured from the stonn top 
downward into the rain. Several features of the plot are wjrth noticing. In the first several range 
bins the rain return power increases because of the progressively greater fractional volume filled 
with predpitaticn. After the peak rai' turn has been reached, corresponding to the first fully 
filled bin, it begiru decreasing with a slope proportional to the attenuation factor. This behavior 
continues until the range bin of the main lobe intercepts the surface; thereafter, the rain return falls 
off more rapidly because of the smaller volume of rain contributing to the backscattered power. 
The surface dutter follows a somewhat reversed progre«ion. This contribution arises from the 
aiuniJus formed by the intersection of con^ntric spheres with a flat earth. As the radar pulse 


10 


propagates through the storm, the annulus expands until a portion of it intersects the main lobe. 
This inteisection corresponds to the broad maximum of the ground return shown in the figurt . 

To describe the sources of error in the rain rate estimates, we first note that in ( 1 7) and (20) 
the measured reflectivity factor at the jth bin, was defined to be 


“ **Rj 

However, the rain return from the jth bin, is not a measured quantity but rather the total 
return power from that bin, We can also account for possible calibration errors by writing 

^R ( ^ ^R ) where is the relative error in and is the tme value. Defining as the 

estimated value of Z^^ ^ then 




( 22 ) 


where Pj is given by (5); i.e.. 


P, = P« + Pf. + Pv, 

Ij Rj Cj Nj 

In computing the statistics of is replaced everywhere in (17) and (20) by 

To characterize the errors in the return power from the precipitation, we write 

Pr = Pr f(K) (23) 

where P|^ is the mean value of P|^ which is given by equation (A4). f(K) is a random variable of 
mean 1 and variance 1/K, where K is the effective number of independent samples [ 18-20] . If we 
stipulate that the radar is to provide contiguous coverage, then for the fairly high spatial resolution 
that will be considered, K is usually too small to guarantee an accurate measurement of Pg . Seveial 
strategies are available, however, to increase K which include: pulse compression [21 ] , the broad- 
band radar ( 22] , frequency agility and multiple beams. An evaluation of these various schemes is 
beyond the scope of this paper. We assume for simplicity that a sufficiently large K is achieved by 
the use of multiple beams. The total peak transmitted power then becomes mPj where m is the 
number of beams and Pj is the peak transmitted power per beam. 


11 


( 24 ) 


In a similar manner, the receiver noise power. . is written 

where K' is the total number of samples per ran^ bin (assuming that all samples of the noise power 
ate independent) and where 

Pj,ak^TBF 

whore is Boltzmann’s constant, T is the receiver temperature in degrees Kelvin, F is the noise 
figure of the receiver, and B is the receiver bandwidth. We have chosen T ■ 290'’K, F ■ S db, and 
B ■ 1/r (where r is the pulse duration) throughout. 

The return powers from the bround, P^^, P^^ are proportional to the backscattering coef* 
ficients Oi", aj". respectively. In general is a function of wavelength, incidence angle, and 
polarization of the incident wave. In the context of the simulation, 0 |°, 02 ° ate chosen to be 
independent and identically distributed log>normal randmn variables. The mean of is taken from 
the vegetation model of Ulaby [23| who h,:s fit the extensive Kansas data to an eight parameter 
equation. The icstUts of the simulation are presented for several values for the standard deviation 
of a°, written S^, which are: • 0.9, 1,8, 2,7 db. The value ■ 1,8 db corresponds approxi* 

mately with that given by Ulaby. For all examples, we have chosen the incident wave to be 
horizontally polarized. 

Substituting (22) and (23) into (2 1) then 

^noj “ *R ^ 


where 


%-PR/(K)tP'^^P,fiK') 


( 26 ) 


The random variables f(K). f(K'), and P^^ are assumed to be mutually independent. 

The estimate of the atteniiation factor is found from (4) where P, is given by (25) with n 
replacing j and where 


To compiete the characterization of the sources in error in y , \ we must 

approximate the errors arising from the use of the k>R, Z-R, and k<^ relationships in these estim- 
ates. Writing, as before, 

k«7R^ 

R-aZ*> 

k«aZ<* 

we make the following assumptions. In the k-R relationship, we fix ( and chose y to be nonnally 
distributed with an unbiased mean and st.<«ndard deviation given by Atlas and Ulbrich ( 14] . 
Similarly, for the Z>R relationship wc fix b and choose a to be nonnally distributed with an 
unbiased mean and a standard deviation equal to 25% of the mean value. The above relationships 
are not independent; therefore (tfi can be expressed in terms of 7 ,|,a,b. i.e., 


d-fb 

(28) 

a* 7,1 

(29) 


Before presenting the results, it is instructive to recognize some of the requirements that , Pq 
must satisfy so that reasonably accurate measurements of fe^y , Rj'^^ can be made. Since 

the most critical variables in the rain rate estimates are and Z^ , we discuss these exclusively. 

To obtain A„ accurately, we must insure that at one or more of the surface bins the following 
inequalities hold: P^ ' » P^ , P(j ' » Pr . If these conditions are satisfied, then from ( 1 ), (2) and 
(4H6) 

An /Pc„ (30) 

where A^ is the true attenuation factor out to the nth range bin. A second necessary condition for 
an accurate estimate of A^ is seen immediately from (30): only if the random variables Uj ° . Oj® 
are dose to one another in a statistical sense will A^ be an indication of the attenuation factor 
rather than the spatial or temporal fluctuations in o'*. Thus, in so far as possible, the incidence 
angle and wavelength should be chosen so that effects such as soil moisture and surface roughness 
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are minimized. This can be partially realized by a selection of K-band frequencies and incident 
angles pointed away from nadir. 

Ak 

From (S). ( 1 3), (22) it follows that an accurate determination of Z requires that P. » 

j j 

Pj, and Pg j “ i. n* This second inequality, however, is in direct opposition to the 

above requirements that P^^ » P|^ at one or more of the surface cells. Since the latter condition 
is essential to the accuracy of the method, the accuracy of must be sacrificed at these surface 
cells. At the precipitation cells, however, the condition P|^ » P[; can be retained. 

It is not difficult to select system parameters so that most of the above conditions are satisfied 
over some band of rain rates. By referring to Fig. 2, it can be seen that P^ » P^ , P^ » P^, over 
all precipitation cells, excepting those at the storm top where a significant portion of the bin is 
unfilled. Furthermore, at several of the surface bins, the conditions Pq » P^ , Pq » P^ hold. 

We conclude that at X = 1.24 cm., R = S mm/hr and for the parameters listed in Table 1, the 
necessary conditions for rain rate estimation via the SRT are fulfilled. 

Results of the simulation are presented in Figs. 3 through 6. Each plot gives the mean and 
standard deviation of R^ y (both normalized by the true rain rate, Rj. ) and of R^ * \ R^ ^ ^ as func* 
tions of distance into the storm for a fixed values of wavelength and rain rate. In all cases, the true 
rain rate is assumed to be constant within the slab extending from z » 0 to z = and zero else* 
where. Since R^ y is not a function of range, we obtain single values for the normalized mean and 
standard deviation. In the figures, this normalized mean is represented by a horizontal line. 

For example, in Fig. 3, where X = 1.24 cm, R » 5 mm/hr and » 1.8 db, the normalized 
mean, R^y/Rp, equals 0.94; that is, the estimate has a small negative bias. At the end of this line 
and corresponding to the nth range bin a vertical bar is plotted of magnitude equal to twice the 
normalized standard deviation of R^ y , 2S^ y , where y = y /Rj . Again referring to Fig. 3, 
we can measure this length by projecting it onto the vertical scale. We fmd that o^ y = 1 .5. The 
same conventions are used to display the normalized mean and standard deviation of R^‘ ^ 

In Fig. 4, for example, at r« 3 km we find that R'* VRx *0.87, o^/Rj ® 0.07 and R'^’/Rj * 0.75, 
Oj/Rp *0.23. 
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As an aside we note that the figures indicate that only in the mid*ranges from about 2,5 to 
6.5 km (comprising about 50% of the total number of range bins) can we even attempt to 
adequately estimate the range-profiled rain rate. This mid-range consists of those bins in the main 
beam that do not intercept the ground and that are also nearly or fully filled with rain. It is easy to 
show that as the antenna diameter to wavelength ratio, D/X, is increased, the number of partially 
filled cells decreases. For example, at DA * 1000, approximately 80% or the bins are relatively 
free from such errors. Such large values of D/X, however, introduce problems of their own. At 
more modest values of D/X, steps might be taken to estimate partial beam-filling (24) or attempt to 
separate, by means of Doppler processing, the surface from the rain return. 

If the wavelength is fixed at X « 1.24 cm and the rain rate is decreased below 5 mmAr, R^y 
soon becomes larger than the tme rain rate (positively biased) while the normalized standard devia- 
tion of R^y becomes progressively larger. For example, at R » 1 mm/hr, R^y ° 1>68 and o^y 
1.1, indicating that the SRT is not feasible at this rain rate and wavelength. By going to rain rates 

s 

larger than 5 mm/hr, the normalized standard deviation of R^ y is significantly reduced while 
accompanied by only a slight increase in negative bias. These trends continue up to about 20-25 
mm/hr where the noise power, Pj^ , becomes comparable to Pq . For rain rates beyond 25 mm/hr, 

n 

the bias in R^y , R^^\ R^^^ begins to increase rapidly. Fig. 4 is plotted for the same parameters 
as in Fig. 3 except that Ry » 20 mm Ar. The increased bias and the reduction of the normalized 
variance with larger rain rate is evident on comparing the two figures. 

If we accept as adequate those estimates with an offset and standard deviation less than 25% of 
the true value, then at X » 1.24 cm, ® 1.8 db, and with the values of Che radar parameters in 
Table 1 , we find that we can estimate the average rain rate between about 7 mm/hr and 20 mm Ar. 

If is reduced to 0.9 db the minimum and maximum rain rates become 4 mm/hr and 25 mm/hr. 
These values and the results at two other frequencies are presented in Table 2. The main determin- 
ing factor for these minimum and maximum rain rates is the magnitude of the total attenuation up 
to the nth range bin. When Ae attenuation, J k(s) ds, is slight (low rain rates) y will 
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generally be greater than the true value: this is due to the logarithm in ( 10 ) which biases values of 
near unity toward the larger values of y . On the other hand, when the attenuation is large 
(high rain rates) the noise power, P» , eventually predominates over . giving from (4H6) an 

n 

attenuation factor A„ larger than the true value. This leads to an underestimate of the attenuation, 
and therefore, an underestimate of the true average rain rate. 

The statistics of \ in the mid-ranges follow the same general behavior as those of 
R^ Y . Notice, however, in Figs. 4 and S that the mean values of R^^ \ R^^ ^ exhibit a gradual 
decrease from about 3 km to 6 km followed by an increase around 6.5 km which corresponds to 
the intersection of the range bin with the surface. Both effects can be explained by the presence of 
the ground clutter which leads to an overestimation of as is shown by (21). This overestimation 
of Z is especially large at the surface cells and implies that for j ^ n, the ratio /S„ is smaller than 
the true value. This, in turn, negatively biases the estimates of (17) and (20). On the other hand, 
as j approaches n, then S^/S^ -*■ 1 and the greatest contributor to the bias is caused by the Z^ ^ 
factor in ( 1 7) and ( 20 ) which produces a positive bias in the rain rate estimates. Tiese effects 
account for the abrupt change in the mean value at about 6 .S km from negatively to positively 
biased values. 

Despite the similarities between R^^ \ R^^\ there are several differences between them. As 
mentioned before, R^* * is independent of errors in the calibration constant . Thus in the usual 
situation where the errors in are 2 db or greater, R^^^ is preferable to R^^^ To show this 
graphically, we plot in Figs. S and 6 the statistics of R^y , R^* \ R^^^ for Rj > 20 mm/hr and 
X " 1.87 cm. In Fig. S the error in Cg is chosen to be zero. In Fig. 6 , C,^ is in error by 3 db; Le. 

C|^ ■ 2 in (2 1). Notice that neither R^ y or R^ ‘ ^ is affected by this offset. 

Although the figures indicate that the variance in R^^ ^ is generally smaller than R*^\ this may 
be somewhat misleading because of the simplified error model used for the k-R, Z-R, k-Z relation- 
ships. In particular, we And from (27), (28) that > » a and 17 ■ /J b'* so that ( 1 7) can be 


written 



J'(l-(l-\^Sj/S,)■‘^'' 
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where 


J'*((l-A„^/O.46 70hS„)*'’’. 

In this form, ^ is no longer explicitly dependent on a and a. Since t}, b, 0 were assumed to be 
Axed, the only sources of variability besides S^, is in y. To obtain a better indication of 
the variance in ^ ^ and a more sophisticated model is required for the characterization of 
errors in the lc*R, Z-R, k-Z relationships. 
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CX)NCLUSIONS 

By using the surface as a reference target, the total attenuation method can be extended to 
spacebome geometries. This method provides estimates of path*averaged rain rate that are inde* 
pendent of the Z*R law and radar calibration constant but dependent on the backscattering coeffi- 
cient of the surface and the k*R relation. Range-profiled estimates of rain rate can be deduced by 
means of a reformulated version of the Hitschfeld-Bordan algorithm. This reformulation was done 
by using the total attenuation, as determined from the suriace reference method, to derive a 
parameter which bounds the error in the original estimate. Of the two estimates that were analyzed, 
the R^ ^ ^ estimate was preferable because of its insensitivity to radar calibration errors. 

To generate the results that have been presented, a simple meteorological model has been used. 
Although the errors caused by simple partial beam-filling effects were included, no spatial variation 
in the rain rate was assumed. Furthermore, the effects of atmospheric absorption, clouds and 
precipitation above the melting layer were neglected as was the presence of non-spherical drops. 
These assumptions were made to decrease computation time and to simplify the presentation of 
the results. Nevertheless, these effects would have to be accounted for in a more rigorous assess- 
ment of the technique. 

In the modelling of the fluctuations in o°, Le. the variability introduced by substituting the 
actual for an equivalent target, we have used in one set of calculations, a » 1 .8 db taken from 
Ulaby ( 23 ] . This value, however, includes data taken over different crop types, soil conditions 
and times of the growing season. Since some of these sources of error will be absent in the data 
used for estimating the attenuation factor, we expect, on the average, to encounter lower values of 
S^. Selected runs through Skylab data [15] taken at X » 2.16 cm for an incident angle of 33*’ have 
yielded a between adjacent 10 km footprints of about 0.43 db. The standard deviation in 
from values separated by S to 10 instantaneous fields of view (IFOV) was found to be 0.5 to 0.6 
db. These values, however, are probably too low for the type of application envisaged since the 
Skylab IFOV (about 4 times the area of that considered in the simulation) would tend to smooth 
over anomalous features appearing in the footprint. Furthcnnore, in the SRT, the changes in soil 
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moisture conditions and the presence or absence of moisture on crops, leaves, etc., would tend to 
make a" more variable than on the average. Because of our lack of knowledge of , the Simula* 
tion was carried out for three values of this quantity. Ttie results have shown that the minimum 
tain rate that can be estimated accurately is a fairly sensitive Ainction of this parameter. 

One of the limitations of the surface reference technique appean to be the narrow range of 
rain rates over which a single frequency radar can work well. In particular, it was found that below 
a certain rain rate, the total attenuation is too small for the technique to work effectively. In this 
region of sli<»ht attenuation, however, the iterative algorithm or the Hitschfeld>Bordan algorithm 
works best. It appears likely that the dynamic range can be extended by using one of the standard 
attenuating wavelength techniques for estimation at the low rain rates. For a dual*wavelength radar, 
additional hybrid techniques would be available. 

As a final comment we remark that the SRT has been analyzed only for rain over land. Over 
ocean it is expected to be less accurate because of the dependence of on wind speed for angles 
greater than 10°. For measurements over ocean however, a rada^radiometer sensor might be 
feasible. In this mode the total attenuation would be obtained from the radiometric intensity 
measurement while tlie radar would provide a measure of the storm height and th<r ran^profiled 
reflectivities. The rain rate algorithms of (1 7) and (20) would then be applicable to this sensor if 
the quantity were interpreted as the radiometrically<derived attenuation factor. 
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APPENDIX 


AVERAGE RADAR RETURN FROM THE RAIN AND GROUND 


In the text, the explicit foims of the mean retum power from the rain and the surface are 
not given. We wish here to write down these expressions and indicate how they .ma;^ be simplified 
uruci »rtain conditions. We begin with an equation for the retum power as a function of time, t, 
from an ensemble of scattereis. A slight generalization of the equation given in (1 ] can be made by 
keeping the attenuation factor i .de the integral, Le., 


p,(t) 




X»P^(t-2r/c) 


- 0.46 

q G^e 



2 

k(i) dt 


dr dS2 


(Al) 


where q is the reflectivity, P^ is the transmitted power, r is the range measured from the radar, X 
is the radar wavelength and G is the antenna gain fiincticn. We choose spherical coordinates (r, 
d, d) with respect to an origin located on the radar at an altitude above the ground. The polar 
angle $ is measured with respect to the z axis which points in the nadir direction. As in the text, 
k(s) is the attenuation coefficient in db/km. In general, q, Sj . S. are functions of (r, 0, d) because 
of the finite extent of the region of precipitation. The antenna gain function G is taken to be 
rotationally symmetric about the direction of maximum gain, d « , d ■ 0, so that, 

G(d,d)-6(d) 

where d is the angle between the vectors (r sin d^, 0, r cos d^) and (r cos d sin d, r sin d sin d, 
r cos d) and is given by 

d » Cos*' (cosd d sin d]^ cos d cos d,g ) 

To simplify (Al) we make the following assumptions: we choose a uniform rain rate from the 
ground up to an altitude z,. We further assume a rectangular transmitted pulse of duration r: 

the retum power is then sampled at times t^ ■ Zr^/c + r, j • 1 n where r.^ , - Tj ■ h ■ cr/2. The 

quantity. Pg (t ■ 2t^Ic+t), is the retum power from the rain at the jth range bin which we denote 

by P. . The final approximatitm is to neglect the rain icFjm frmn all but the main lobe. Equation 

j 

(Al) becomes 
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fr^‘l f G. *(♦■«) F(r,».«) ^ ^ 


J (4ir) 
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t| 4« 


dR 


(A2) 


whw 

G wirhin main lobe 

c«- 

0 ebewheie 

and F is a nn factor, being unity for that region of the range bin which intercepts the rain and zero 
elsewhere. For a fixed $, s^ • (z^ • z,) sec 0 , Sj * r so that s, • S| is simply the distence through 
the stonn to range r at a particular value of d. Using the relation f) ■ IK 1^ ZfK* 1 2] . then 


f'j»h 


«CZ f f G„^(d,d)F(rAd) e 
i ^ 


-046 (r - Attcd) 


dr dl2 


(A3) 


4X t; 


Where, 


C*»^ IK|2 P^/64X* 

A ■ (Zo - z,) 

lKl^-jm2 4-2|/|m*-ll 

m • complex index of refraction of water 

The r integration in (A3) can be done for an arbitrary range bin; this results, however, in tliree 

different expressions for P. depending on whether the bin is partially Ailed at the storm top or at 

j 

the ground cell of whether the cell is fully Ailed. To be brief, we write down P. only for a bin 

J 

Ailed with rain. The remaining cas« can be done in a similar manner except that greater care is 
needed in i^cifying the d limits of integration. 

For a fully Ailed bin, F * 1. Integrating by parts, then 



— dr • I 

r2 (r^+h) 


•♦•7(Ej(7(r^4>h))- Ejtytj)) 
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where, 


lij(x)= j -p ill 


7 * 0.46 0! 

Assuming that 7 Fj » 1, we can use a large argument expansion on Ej(x) [3] . Omitting «9(rj'^) 


terms then 

e”^ 1 -71. -7h ^ t 

f dt * — e^ J (1 “ e^ rj^/(rj + h)^) 

.j 

The return power from a uniformly filled range bin then becomes 


‘i 7r/ \ L 


(A4) 


If the attenuation per range bin is small, i.e. 7 h « 1, and Tj » h then 

7-*(l-e-T^"(rj/(rj+h))2)5:h 

Furthermore, if the main beam is sufficiently narrow so that sec 6 where is the angle 

between nadir and the direction of maximum gain, then 


p ^-7(rj-Asecg^) 


(A5) 


where 


having used the fact that 


h |KI= Pt ^ 

C = — — — I C-(^) sin di^ 
32 „ 


0b/2 


(A6) 


j G^^(0,d)df2 = 2jrJ S^(0)sin0d0 

4ir 0 

where dg is the angular extent of the main lobe. Thus, (A6), for a symmetric antenna pattern, 
involves only a single integration over 0. 
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F<m‘ the purpose of the simulation in the text, the mean return power has been generated 
using (A4) (and the corresponding expressions at the storm top and bottom where the beam is 
partially filled) instead of (AS) since the assumptions 7 h « 1 , sec 9 ^ sec 9^ are not always valid 
for the frequencies and beamwidths that we wish to consider. 

For a nonuniform rain rate, the r integration can not be done so that a volume integration 
is necessary for finding P. . If the main beam is sufficiently narrow, however, and yh « 1, then 
we have, approximately, for a filled range bin 




(A7) 


where Aj is the attenuation factor ^ven by. 


H).46 h 
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i«l 




(A8) 


and C,j is given by (A6). Defining the measured reflectivity factor by Aj then (A7) 

becomes 






(A9) 


To obtain the backscattered power from the surface, we begin with the expression ( 1 ), 




(4ir)^ 
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P^(t - 2r/c) 0 ° (4,9) 


-0.46 
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k(s) di 


dS' 


(A 10) 


where the region of integration is over the entire ground plane. We again assume a rectangular 
pulse of duration r and sample the return waveform at times tj ® 2r./c + r, j = 1, .... n. The notation 
Pq (tj = 2tj/c + r) “ pQ^ is used. Assuming that h « r^ 9^, where 9g is the angular extent of the 
main lobe, then 9 varies only slightly over the region of integration. This implies that the attenua- 
tion is essentially constant over the area of integration so that it can be taken outside the integral. 

If we furthermore assume that a° is independent of 0 then o®(0,9) ~a®(9^ ), where 9^ is defined 
below. Under these assumptions, the return power from the surface, in the presence of 
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precipitation, is 


''i (4ir)5 

where is given by (A8) and 


Pt- I / — ; TT P 

0«o p»Pj, ) 




Pb»((rj + h)2 - 

Performing the p integration and using the approximation 

1 1 


Tj (tj + h)- 




which is suitable for r^ » h, then 


where 


PG' = Cca“(e„)A/r/ 

=Tan'^((p, +p^)/2Zq) 

2V 

Cc =(X^ ?X h/(4t7)3) / G2(0„,0)d0 


(All) 

(A 1 2) 
(A 13) 


The calculation of the return power from the surface in the absence of precipitation, , follows 
at once from (A1 1) by setting A^ equal to unity, i.e., 

<."(«„ )/r,> 

We now wish to indicate how can be found in terms of measurable quantities. From (A9), 


Zm, ■ ')' P./Cr ‘AH) 

where Pp is the rain return power from the jth range bin. Implicit in the choice of calibration 

‘‘j 

constant C,^ is the assumption of a fully filled range bin. If the extent of the beam filling were 
known, , could be replaced by a calibration constant which takes this into account. However, 


we assume that no such apriori knowledge is available. Thus, the offset errors of in (A 14) 

arise from partial filling of the range bin, the normal calibration errors in , and the additional 

errors introduced by the approximations that led to (A9). Noticing, moreover, that P, rather tlian 

‘j 

P,^ is a measured quantity, (A14) must be changed to 




(A15) 


where from equation (S) of the text, 


» -P. +Pg +Pm 


(A16) 


The terms of P, , other than P_ , are seen to introduce an additional offset into . For the 
‘i "‘j 

p«,»*s of stoul.dc. i. U.. d.. Z„ . J = 1 „ d... .pp« 1. d.0 to. to. 

estimates are generated by means of (AIS) and (A16). The mean value of PL is found from (All) 
through (A 13); the mean of P. , for a range bin uniformly filled with rain, is found from (A4). 
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Table 1 


Pt 

a 

peak transmit power = 1 kw 

h 

a 

range resolution * 0.3 km 

D/X 

a 

antenna diameter to wavelength ratio « 300 

G(0) 

» 

antenna gain = 1 1 



/ 1 Ej ((?) 1^ d n 

Ej (fl) 

= 

incident electric field « Jj (»r D sin (9/ X) /(D sin fl / X) 

g(P) 

s 

antenna aperture illumination » 1 (o < p < 0 / 2) 

PRF 

= 

pulse repitition frequency = 10 kHz 

m 

= 

number of multiple beams 4 


c 

storm height = 5 km 


= 

satellite altitude » 700 km 


= 

satellite velocity = 7 km/sec 
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Table 2 


Wavelength (cm) 

S„ (db) 

Min R (mm/hr) 

Max R (mm/hr) 

0.86 

0 

<1 

12 

f 

0.9 

2 

1.8 

3.5 

2.7 

5 

1.24 

0 

<1 

25 

0.9 

4 

1.8 

6 

2.7 

8 

1 

1.87 

0 

1 

55 

0.9 

8 

1.8 

14 

2.7 

20 
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Figure 1 . Spacebome Radar Geometry 
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